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Abstract 

We calculate the electromagnetic form factors of a bound proton. The Chiral Quark-Soliton 
model provides the quark and antiquark substructure of the proton, which is embedded in nu- 
clear matter. This procedure yields significant modifications of the form factors in the nuclear 
environment. The sea quarks are almost completely unaffected, and serve to mitigate the valence 
quark effect. In particular, the ratio of the isoscalar electric to the isovector magnetic form factor 
decreases by 20% at = 1 GeV^ at nuclear density, and we do not see a strong enhancement of 
the magnetic moment. 



I. INTRODUCTION 

Recent polarization transfer experiments at TJNAF Q observed a difference in the elec- 
tromagnetic form factors of a proton bound in a Helium nucleus compared to a free one. 
This, along with other effects, such as the nuclear EMC effect j^, seems to suggest the 
modification of hadrons in the nuclear medium. 

There is extensive work on the medium modifications of electromagnetic properties of 
the nucleon in the literature (for example, see Refs. This includes effective 

Lagrangians as well as models that include the quark substructure of hadrons. While in 
principle these effects could be couched in terms of effective field theory operators, it is our 
thesis that such results may be more transparent, physically intuitive or straightforward to 
calculate when viewed as a change in the internal structure of the hadrons. 

we wn. use the Ch.a.-Qua.. Sol.ton n,ode. (CQS) flQ. whieh has a direct eon.eetion 
to QCD via the Instanton Liquid model, to provide our subnuclear degrees of freedom. The 
primary motivation is that this model includes sea quarks which we have seen to be im- 
portant in the nuclear EMC effect In that case, the large medium modification in the 
valence quark sector is reduced through the lack of such an effect in the sea (which can be 
seen directly in Drell-Yan experiments The CQS is combined with the nuclear medium 

in a self-consistent quark-meson coupling calculation as in our previous work [9], and the 
electromagnetic form factors are extracted via the wave functions of the quarks using the 
results of Ref. The overall procedure is similar to the Quark-Meson Coupling model 

(QMC) which uses the MIT bag model for the nucleon. The bag model does not include 
sea quarks. It is a confining model, whereas the CQS model is not. Additionally, the QMC 
model calculation, when coupled with a Relativistic Distorted Wave Impulse Approxima- 
tion (RDWIA) calculation \12^ or a Relativistic Multiple-Scattering Glauber Approximation 
RMSGA) calculation [13|, |l4|, improves the agreement between theory the TJNAF data 
^. With our study, we hope to reinforce the interpretation of the medium effect in terms 
of quark degrees of freedom, as well as provide an alternate model when the accuracy of the 
data is improved. 

We begin with a brief description of the CQS model in Section |n] In Section IIIIl we 
motivate and present our procedure to embed this model in nuclear matter. This description 
differs only slightly from that in our previous work 9]; it is repeated for completeness. 
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Subsequently, we describe the numerical methods, and proceed to the results in Section IIVI 



II. CHIRAL QUARK-SOLITON MODEL 

The CQS model Lagrangian with (anti)quark fields ^, and profile function 6(r) is [3] 

£ = ^(i^ - Me*''^^-^®("))^, (2.1) 

where 0(r ^ oo) = and 0(0) = — vr to produce a soliton with unit winding number. 
The quark spectrum consists of a single bound state and a filled negative energy Dirac 
continuum; the vacuum is the filled negative continuum with = 0. The wave functions in 
this spectrum provide the input for the electromagnetic form factors. 

We work to leading order in the number of colors {Nc = 3), with Nj = 2, and in the 
chiral limit. While the former characterizes the primary source of theoretical error, one 
could systematically expand in Nc to calculate corrections. We take the constituent quark 
mass to be M = 420 MeV, which reproduces, for example, the A^-A mass splitting at higher 
order in the Nc expansion, as well as many 

The theory contains divergences that must be regulated. We use a single Pauli-Villars 
subtraction. The Pauli-Villars mass is determined by reproducing the measured value of 
the pion decay constant, = 93 MeV, with the relevant divergent loop integral regularized 
using Mpv ~ 580 MeV. 

The nucleon mass is given by a sum of the energy of a single valence level [E^), and the 
regulated energy of the soliton {Eq, equal to the sum of energy levels. En, in the negative 
Dirac continuum with the sum of the energy levels in the vacuum, En^ , subtracted) 
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= NcE- + EeiM) - -^Ee{Mpv) (2.2a) 



M'py 



(2.2b) 

M=M' 



The field equation for the profile function, which follows from the Lagrangian (j2.H) . is 

e{r) = aictan^-^^, (2.3) 



where and p^^ are the quark scalar and pseudoscalar densities, respectively, and are given 
by sums of the wave functions of every occupied energy level. 



The electromagnetic form factors are also given in terms of the wave functions, and 
are derived in Ref. The formulae are reproduced here, with a Pauli-Villars regulator, 
for convenience. To leading order in Nc-, we have only the isoscalar electric and isovector 
magnetic form factors (G^"'''^ = ± G\) 

Gl-\q') ^-=°° ^ / dre^ €{r)Mr) - ^^^\r)i^^\r)\ (2.4a) 

"J n \ J \ T? ^T7V 



E^<E^ 

E i^r'\r)^Hr'^r\r\ (2.4b) 

The^rV) 

are the solutions of the Dirac equation with the replacement M — ^ Mpy. In the 
nuclear medium, Eqs. ()2.4|) acquire a dependence on the Fermi momentum G^~'''^(g^) — *■ 
fci?) through the wave functions. This dependence is the subject of the next 

section. 

III. NUCLEAR PHYSICS 

We will begin with some motivation for our procedure to couple the quark substructure 
of the nucleon to the nuclear medium. Through the use of QCD sum rules, loffe derived 
a relationship between the vacuum scalar condensate, (V'^)o) and the nucleon mass. One 
can re-derive this estimate in a constituent quark field theory such as we are using here. We 
begin with the scalar condensate 



(^^)o = -tr 



(27r)4;j/ - M 
NcMK^ 



(3.1) 



47r2 ' 

where the divergent integral is regulated by a momentum cutoff (playing the role of the 
Borel mass in the QCD sum rule approach). Using the fact that constituent quarks are 
essentially defined as having a mass ~ M^/Nq, we can rewrite Eq. p.ip as 

47r2 _ 

Mn ~ -j^{iJtP)o. (3.2) 



Although Eq. ()3.2j) is not a very accurate estimate, it does highlight the role of the conden- 
sate. It will be modified in the presence of other nucleons. 

The condensate at finite density can be written in terms of the nuclear scalar density 
and the nucleon sigma term [l6| as 



(^^), = (^^)o - (^^)o^pf . (3.3) 



We can then substitute Eq. (j3.3|) into Eq. (j3.2p to obtain a schematic picture of the effect 
of the nuclear medium on the nucleon mass 

A-k"^ — 

MM ~ [(^V^)o - c.pf ] , (3.4) 

where Cg is the combination of the vacuum condensate, pion mass, decay constant and the 
the sigma term in Eq. ()3.3p . 

Using this dependence of the nucleon mass on the nuclear medium as a guide, we incor- 
porate the medium dependence in the model by simply letting the quark scalar density in 
the field equation ()2.3|1 contain a (constant) contribution arising from other nucleons present 
in symmetric nuclear matter. This models a scalar interaction via the exchange of multiple 
pairs of pions between nucleons. We take the scalar density to consist of three terms: 1) the 
constant condensate value {ipip)o (in the vacuum or at large distances from a free nucleon), 
2) the valence contribution and 3) the contribution from the medium which takes the 
form of the convolution of the nucleon and valence quark scalar densities as in the QMC 
model |3|] 



plir) ^ (V^^)o + p:{r) + csj dr'p^ir - r')pl{r') (3.5a) 
= {i^^), + pl{r) + ~CspfS (3.5b) 
[ dr'plir'). (3.5c) 



We take the pseudoscalar density to have only the valence term p^^ ~ p^^; the two other 
contributions analogous to the first and third terms of Eq. ()3.5|) vanish due to symmetries 
of the QCD vacuum and nuclear matter. These approximations to the densities neglect 
the precise form of the negative continuum wave functions in Eq. ()2.3|1 . The resulting free 
nucleon profile function has no discernible difference from a fully self-consistent treatment, 
demonstrating the excellence of this approximation. We take Cg = Cg/S in Eqs. fj3.4p and 
()3.5|) to be a free parameter, which we vary to fit nuclear binding. This can be seen as either 



varying ctat in Eq. (|3.3p or the vacuum value of the condensate in Eq. ()3.5|) with pp Pr/cs, 
as was done in Ref. j^l, since the overall normalization cancels in Eq. fl2.3|) . 

The nucleon scalar density is determined by solving the nuclear self-consistency equation 

The dependence of the nucleon mass, and any other properties calculable in the model, 
on the Fermi momentum kp enters through Eq. ()3.6|) . Thus there are two coupled self- 
consistency equations: one for the profile, Eq. ()2.3|) . and one for the density, Eq. ()3.6p . 
These are iterated until the change in the nucleon mass Eq. (|2.2j) is as small as desired (in 
our case, AMjy ^ 0.1 MeV) for each value of the Fermi momentum. We use the Kahana- 
Ripka (KR) basis with momentum cutoff A and box size L ex trap olated to infinity 
(from a maximum value of AL = 150, comparable to that in Ref. |ll|), to evaluate the 
energy eigenvalues and wave functions used as input for the densities, nucleon mass, and 
electromagnetic form factors. 

While the vacuum value of the condensate does not vary with the Fermi momentum by 
definition, the effective condensate, (V'V')o + Cgp^ {kF)S{kF), falls ~ 30% at nuclear density, 
q.v. Eq. ()3.3|) . This is consistent with the model independent result |3] that predicts a value 
25-50% below the vacuum value. 

A phenomenological vector meson (mass = 770 MeV) exchanged between nucleons 
(but not quarks in the same nucleon), is introduced as a substitute for uncalculated soliton- 
soliton interactions in order to obtain the necessary short distance repulsion which stabilizes 
the nucleus. This does not affect the form factors Eqs. ()2.4aj) and ()2.4bjl . The resulting 
energy per nucleon is 

T = / 7^VA:2 + M^(M^ + i^4pb(M- (3.7) 

A pBikp) J (27r)'^ 2mi 

The mass of a free nucleon is computed to be Mjq{kp = 0) = 1209 MeV. The ~ 30% 

difference is as expected in the model at leading order in Nc- We evaluate the nucleon 

mass Eq. (|2.2|) and energy per nucleon Eq. (|3.7p as a function oi kp. We choose our free 

parameters to fit E/A — Mjv(O) = B = —15.75 MeV at the minimum. We use the value 

Cs = 1.27 (corresponding to ctat = 41.4 MeV), and vector coupling gl/'ir^ = 10.55, which 

gives a Fermi momentum of kp = 1.38 fm~^ in nuclear matter consistent with the known 



value kp = 1.35 ± 0.05 fm"^ 



We plot the binding energy per nucleon using Eq. ()3.7|l 



in Fig. The compressibility is = 348.5 MeV which is above the experimental value 
K = 210± 30 MeV, but well below the Walecka model Q value of 560 MeV. The self- 
consistent calculation results in the profile functions for zero density, 0.5po, 1-Opo and 1.5po 
in Fig. 121 (where po is nuclear density). 

IV. RESULTS AND DISCUSSION 

We use Eqs. fl2.4aj] and (j2.4bjl to calculate the form factors, which we present in Figs. El 
and 13 We also present the results in terms of the ratios 



These ratios are plotted in Figs. El and El for 0.5po) 1-Opo and 1.5po. 

The electric form factor is dominated by the valence contribution and shows a dramatic 
effect, while the magnetic form factor has equally important contributions from the valence 
and the sea. The latter shows almost no change in nuclear matter; it shows only a 1.3% 
enhancement of the magnetic moment at nuclear density, and a 2.3% enhancement at 1.5 
times nuclear density. The effect in the electric form factor calculated here is comparable 
to that of the QMC model '3]; the main difference from that calculation lies in the lack of 
enhancement in the magnetic form factor, specifically the practically unchanged value of the 
magnetic moment. 

While both form factors use the same wave functions, the isovector magnetic form factor 
includes an extra weighting by a factor of the angular momentum of the state (relative to 
the electric form factor) due to the in Eq. ()2.4b|) . This extra factor is not only responsible 
for making the regularization of Eq. ()2.4b|l necessary, but for making the sea contribution 
as important as the valence. In the CQS model, the orbital angular momentum carried by 
the sea is comparable to the orbital angular momentum carried by the valence quarks [20| 
(the sum of which make up about 60% of the total angular momentum of the nucleon state, 
with the remainder belonging to the intrinsic spin of the constituent quarks) . 




(4.1) 



where -g^ = Q2^ X is E{T = 0) or M(T = 1), and the double ratio 

GUQ')/G*m{Q') 

GEm/GMm' 



(4.2) 
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FIG. 1: Binding energy per nucleon B = E/A — Mi\j. The box and the gray band correspond to 
the uncertainty in the known values of the binding energy, density and compressibihty of nuclear 
matter 
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FIG. 2: Profile functions in nuclear matter. The solid line is the profile function for 1.5poi the 
curves with progressively longer dashes correspond to l.Opoi O.Spo and zero density respectively. 
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FIG. 3: The isoscalar electric form factor at nuclear density (solid) and at zero density (dashes). 
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FIG. 4: The isovector magnetic form factor at nuclear density (solid) and at zero density (dashes). 

Conversely, the isoscalar electric form factor (which is finite, after the vacuum subtrac- 
tion) does not have as large of a contribution from the sea. The valence level is the most 
important piece, even at > 0, since the dependence in the form factors arises from 
the wave functions The negative Dirac continuum wave functions largely cancel in the 
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FIG. 5: The electric (lower three curves) and magnetic (upper three curves) form factor ratios in 
Eq. (|4.1() for O.bpo (long dashes), l.Opo (solid) and 1.5po (short dashes). 
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FIG. 6: The double ratio Eq. H4.2|) of the electric to magnetic form factors in nuclear matter and 
in the vacuum from the CQS model (heavy) and the QMC model 3] (light). Three densities are 
shown: O.Spo (long dashes), l.Opo (solid) and 1.5po (short dashes). 
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vacuum subtraction in Eq. ()2.4a|) . 

The magnetic form factors are sensitive to the tail of the quark wave functions, and 
the mere existence of a tail is due to the lack of confinement. This is one reason for the 
discrepancy between the current results and the QMC model P], but the primary source 
is due to the resistance to change of the sea. The former accounts for only a few percent 
of the difference; it is the latter that is our most important result. We see that the role of 
antiquarks is again prevalent as in our previous work j^. 

The double ratio obtained in Fig. IHlhas the same trend as the QMC model 0], but differs 
in the details. Since we obtain a similar double ratio, we expect to have similar results if we 
compare these results with the polarization transfer data This requires one to take the 
final state and relativistic effects into account through the use of the RDWIA lial or the 

n 

RMSGA [1J|, which accounts for a few percent of the discrepancy between the results for 
bound and free protons. A RMSGA calculation for the Helium reaction studied in Ref. 

n 

has been done with these CQS model results |21|], and it delivers remarkably similar results 
to the same calculation done with the QMC model jlj]. The CQS model predicts a smaller 
deviation than the QMC model from a Relativistic Plane Wave Impulse Approximation 
(RPWIA) calculation, which is taken as a baseline in Ref. While it slightly worsens the 
agreement with the data at ~ 1, the differences are of the same order of magnitude as 
the current experimental error, and both models under predict the observed deviation from 
a RPWIA calculation. At higher Q^, the two models produce nearly identical results for 
Helium. 

We ignore important corrections due to the rotation of the soliton that are suppressed 
by 1/Nc- These corrections break the N — A degeneracy, and improve the agreement of 
the vacuum form factors with experiment [lli. More relevant to the calculation presented 
here, these corrections do not affect the dependence, but instead affect the normalization 
of the form factors [11]. However, there is no reason at that level to continue to ignore 
quantum fluctuations of the the pion fleld (quark loops), and treat the proflle function as a 
purely self-consistent mean fleld. We will save this difficult problem for the future. 

We have calculated the electric and magnetic form factors at leading order in Nc at 
nuclear density using the CQS model. Our results help validate the apparent success of 
the QMC model in describing the polarization transfer experiment 1^, and provide a 
counterpoint to be distinguished when finer resolution becomes available in the data. In 
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fact, the difference between the CQS model double ratio and the QMC model j3| is roughly 
the size as the current experimental error. Specifically, data on the bound nucleon magnetic 
form factor at low Q^, particularly the magnetic moment, could serve to determine the role 
of sea quarks in nuclei. 
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